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Abstract 

A generalization of the Apery- like numbers, which is used to describe the special values Cq(2) and C<p(3) of 
the spectral zeta function for the non-commutative harmonic oscillator, are introduced and studied. In fact, 
we give a recurrence relation for them, which shows a ladder structure among them. Further, we consider 
the 'rational part' of the higher Apery-like numbers. We discuss several kinds of congruence relations among 
them, which are regarded as an analogue of the ones among Apery numbers. 
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(N ■ 1 Introduction 

> 

The non- commutative harmonic oscillator is the system of differential equations defined by the operator 
where a and (3 are real parameters. In this paper, we always assume that a > 0, (3 > and af3 > 1. Under these 



o 



conditions, one can show that the operator Q defines an unbounded, positive, self-adjoint operator on the space 
L 2 (M;C 2 ) of C 2 -valued square integrable functions which has only a discrete spectrum, and the multiplicities 
m(A) of the eigenvalues A £ Spec(Q) are uniformly bounded [22]. Hence, in this case, it is meaningful to define 
its spectral zeta function (q(s) — TrQ~ s — J2\eSpcc(Q) m (A)A~ s . This series converges absolutely if 5Rs > 1, 
and hence defines a holomorphic function on the half plane -Jis > 1. Further, Cq( s ) is meromorphically continued 
to the whole complex plane C which has 'trivial zeros' at s = 0, —2, —4, . . . (see [7], [21]). 
The aim of this paper is to study the higher Apery-like numbers Jk (n) defined by 



J fe (n) := T 



(1 — x\ )(1 — x\ ■ ■ ■ xt) \ n dx\dxi ■ ■ ■ dxk 



[(Ml* V (i - %i ■ ■ ■ H) J 1 



for k > 2 and n > 0, which are a generalization of the Apery- like numbers J2(n) and Js(n) studied in [11]. This 
object arises from the special values of the spectral zeta function (q(s): In [8], the generating functions of the 
numbers foin) and Jz{n) are used to describe the special values Cq(2) and Cq(3) of the spectral zeta function 
Cq(s). Similarly, the higher Apery- like numbers Jfe(n) are closely related to the special values Cq(^) ( see §3-3). 

We first show that Jk{n) satisfy three-term (inhomogeneous) recurrence relations, which is translated to 
(inhomogeneous) singly confluent Heun differential equations for their generating functions. The point is that 
these relations or differential equations are connecting Jfe(n)'s and Jfc„2(^)'s. This fact implies that there could 
be a certain relation between Cq(^) an d Cq(^ ~ 2). It would be very interesting if one can utilize these relations 
to understand a modular interpretation of Cq(4),Cq(6), ■ • ■ based on that of (q(2) (see [12]). We also notice 
that these recurrence relations quite resemble to those for Apery numbers used to prove the irrationality of C(2) 
and C(3) (see [27]), and this is why we call Jfe(n) the (higher) Apery-like numbers. 
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By a suitable change of variable in the differential equation, we also obtain another kind of recurrence 
relations, which allow us to define the rational part of the higher Apery-like numbers (or normalized higher 
Apery-like numbers) Jfc(n). In fact, each Jfe(n) is a linear combination of the Riemann zeta values ((h), — 
2), . . . and the coefficients are given by J TO (n)'s. Since there are various kind of congruence relations satisfied 
by Apery numbers (see, e.g. [4], [5], [1]), it would be natural and interesting to find an analogue for our higher 
Apery-like numbers. Actually, we give several congruence relations among J 2 (n) and J^{n) in [12]. We add 
such congruence relations among Jk(n), and give some conjectural congruences. 

2 Apery numbers for C(2) and ((3) 

As a quick reference for the readers, we recall the definitions and several properties on the original Apery 
numbers. 

2.1 Apery numbers for £(2) 
Apery numbers for £(2) are given by 

fc=0 v 7 v 7 k=0 v 7 v 7 V m=l m=l W \ m J / 

These numbers satisfy a recurrence relation of the same form 

n 2 u{n) - (lln 2 - lln + 3)u(n- 1) - (n - l) 2 u(n - 2) = (n > 2) (2.1) 

with initial conditions A 2 (0) = 1,^(1) = 3 and B 2 (Q) = 0,^(1) = 5. The ratio B 2 (n)/A 2 (n) converges to 
C(2), and this convergence is rapid enough to prove the irrationality of ((2). Consider the generating functions 

oo oo 

A 2 (t) =J2Mn)t n , B 2 (t) = J2 B 2(n)t n , TZ 2 (t) = A 2 (t)((2) - B 2 (t). 

n=0 n=0 

It is proved that 



L 2 A 2 {t) = 0, L 2 B 2 (t) = -5, L 2 U 2 {t) = 5, 



where L 2 is a differential operator given by 

d 2 d 
L 2 = t(t 2 + lit - 1)— + (3t 2 + 22t - 1)— + {t + 3). 

The function lZ 2 (t) is also expressed as follows: 

dxdy 



n 2 {t) 



JO 



xy + txy(l -x)(l - j/)' 



The family Q\ : 1 — xy + txy(l — x)(l — y) = of algebraic curves, which comes from the denominator of 
the integrand, is birationally equivalent to the universal family C 2 of elliptic curves having rational 5-torsion. 
Moreover, the differential equation L 2 A 2 (t) = is regarded as a Picard-Fuchs equation for this family, and 
A 2 {t) is interpreted as a period of C 2 (see [2]). 
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2.2 Apery numbers for £(3) 
Apery numbers for £(3) are given by 

2/ , t \ 2 ra / x 2 / , , \ 2 



i- — n \ / \ / h—(\ \ / \ / \rn = l m = 



' 2m 3 ( ™ )( n+m ) ] 

fe=0 x 7 x 7 fe=0 x 7 x 7 Vm=l m=l Vm/ V m I / 

These numbers satisfy a recurrence relation of the same form 

n 3 u(n) - (34n 3 - 51n 2 + 27n - h)u{n - 1) + (n - l) 3 w(n - 2) = (n > 2) 

with initial conditions ^3(0) = 1,^4 3 (1) = 5 and B 3 (0) — 0, ^(l) = 6. The ratio B 3 (n)/A 3 (n) converges to 
C(3) rapidly enough to allow us to prove the irrationality of ((3). Consider the generating functions 

00 00 
A s (t) = J2Mn)t n , S 3 (i) = ^S 3 (n)i", Tl 3 {t) = A 3 (t)((3) - B 3 (t). 

n=0 n=0 

It is proved that 

L 3 A 3 (t) = 0, L 3 B 3 (t) = 5, L 3 U 3 (t) = -5, 
where L 3 is a differential operator given by 

d 3 d 2 d 

L 3 = t 2 (t 2 - Mt 2 + 1)— + t(6t 2 - 153* + 3)^ + (7t 2 - 112* + 1)— + (* - 5). 



The function 1Z 3 (t) is also expressed as follows: 



Jo Jo Jo 1 



dxdydz 



(1 - xy)z - txyz{\ - x)(l - y)(l - z) ' 



The family Q 3 : 1 — (1 — xy)z — txyz(l — x)(l — y)(l — z) = of algebraic surfaces coming from the denomi- 
nator of the integrand is birationally equivalent to a certain family Cf of K3 surfaces with Picard number 19. 
Furthermore, the differential equation L 3 A 3 {t) — is regarded as a Picard- Fuchs equation for this family and 
A 3 {t) is interpreted as a period of C 3 (see [3]). 

2.3 Congruence relations for Apery numbers 

Apery numbers Ai (n) and A 3 (n) have various kind of congruence properties. Here we pick up several of them, 
for which we will discuss an Apery-like analogue later. 

Proposition 2.1. Letp be a prime andn = n +nip+- ■ -+rikp k be thep-ary expansion ofn € Z> (0 < nj < p). 
Then it holds that 

k k 
A 2 (n) = Y[A 2 (n 3 ) (mod p), A 3 (n) = J| A 3 {nj ) (mod p). 

3=0 j=o 

Proposition 2.2 ([4, Theorems 1 and 2]). For all odd prime p, it holds that 

A 2 (mp r - 1) = A 3 (mp r ^ 1 - 1) (mod p r ), 
A 3 (mp r - 1) = A 3 (mp r ~ 1 - 1) (mod p r ) 

for any m,r G Z >0 . These congruence relations hold modulo p 3r if p > 5 (known and referred to as a super- 
congruence) . 
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We denote by tj(t) the Dedekind eta function 

oc 

V(r)=q 1/24 l[(l-q n ), 9 = ^" (9r > 0). (2.2) 



n=l 



Proposition 2.3 ([28, Theorem 13.1]). For any odd prime p and any m,r € Z>o urai/i m orfrf, it Zio/rfs i/iai 

A 2 (H^i) - A p A 2 ( mp " 2 1 - 1 ) + (-1)^- 1 )/V^ 2 ( mp " 2 2 ' 1 ) = (mod/). (2.3) 
/fere A„ is defined by 

f^X n q n =V^rf=qf[(l-q 4n f. 

n—1 n—1 

Proposition 2.4 ([5, Theorem 4]). For any odd prime p and any m,r G Z >0 itf&ft m & ftoZds £fta£ 

^3( ! ^)-7 P ^ 3 ( !:i£: ^)+P 3 A3(^ 2 ^)^0 (mod/). (2.4) 

Here j n is defined by 



J2 lnq n = ry(2r) 4 ?? (4r) 4 = q J] (1 - g 2 ") 4 (l - a 4 



n—1 n—1 



3 Apery-like numbers for Cq(2) and (q (3) 

We introduce the Apery like numbers J2(n) and J 3 {ri), and give a brief explanation on their basic properties 
and the connection between the special values Cq(2), Cq(3) of the spectral zeta function (q(s). 

3.1 Definition 

We define the Apery-like numbers for (q(2) and (q(3) by 

rl r 1 ' {I - x\){l - xi)\ n d Xl dx 2 



J 2 (n) :=4 



./ n \ (1 j 2 / 1 £l#2 



Js(n) ;=8 ( {1 -4)(l-44) \ n d Xl dx 2 dx 3 

JO JO Jo V (1 — X 1 X 2 X 3) 2 / 1 — x \ x \ x \ 

The sequences {J2M} and {J 3 (n)} satisfy the recurrence formula (Propositions 4.11 and 6.4 in [8]) 

An 2 J 2 {n) - (8n 2 - 8n + 3) J 2 (n - 1) + 4(n - l) 2 J 2 (n - 2) = 0, (3.1) 

2 n (n — IV 

4n 2 J 3 (n) - (8n 2 - 8n + 3) J 3 (n - 1) + 4(n - l) 2 J 3 (n - 2) = — ^ -f- (3.2) 



(2n- 1)!! 



with the initial conditions 



J 2 (0) = 3C(2), J 2 (l) = |c(2); J 3 (0) = 7C(3), J 3 (i) = ^C( 3 ) + ±. 

It is notable that the left-hand sides of these relations have the same shape. Since the relations (3.1), (3. 2) and 
the one (2.1) for A 2 (n) have quite close shapes, we call the numbers J 2 (n) and J 3 (n) the Apery-like numbers. 
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3.2 Generating functions and their differential equations 

The generating functions for J2(n) and Jz{n) are defined by 

00 /■! yl ^22 

«*(*) := = 4^ y (^^J^ 1 ! 2 ^,^ ^, (3-3) 

~ f 1 f 1 f 1 1 - x 2 x 2 x 2 

By the recurrence relations (3.1) and (3.2), we get the differential equations 

V m w 2 {t) = 0, (3.5) 
V H w 3 (t) = \ 2 FA 1,1; ] , (3.6) 



2" * V 2 

where £> H denotes the singly confluent Heun differential operator given by 



(3.5) is solved in [18] as 



A, = t(l - tf^ + (1 - 3f)(l - i)| + < - 2. (3.7) 



1 -t V 2 2 *- !, 
Here 2^1 (d, 6; c; z) is the Gaussian hypergeometric function. Now it is immediate that 

'n 



1=0 

Similarly, (3.6) is solved in [11] as 



J 



Mn) =3C(2) 3 . . (3.8) 



00 / n / ^ ^ fc— 1 -j / 1\ — 2 

2 



Therefore it follows that 



W2i + l) 3 V 7 

n=0 \fe=0 x / \ / j=Q v J ' \ ■' 



1\ 2 / \ " / 1\ 2 / \ J" 1 1 / l\-2 



Remark 3.1. The function 



^2(^=2^1 I,I;1;T 2> ) =-^(1-*)^) (T 2 



2'2' ' y 3C(2) V ; zv; V *-l 



satisfies the differential equation 



T(T ' " + (3T2 " + T ) W2{T) = °' 

which can be regarded as a Picard-Fuchs equation for the universal family of elliptic curves having rational 
4-torsion [12]. This is an analogue of the result [2] for the Apery numbers for £(2) (see also Section 2.1). It is 
natural to ask whether there is such a modular interpretation for Ws(t) (or "W^T)"). We have not obtained 
an answer to this question so far. ^> 
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3.3 Connection to the special values of (q(s) 

We also introduce another kind of generating functions for Jfe(n) as 



92( 



dxi dx< 



oio V /(l-.Tf a; 2)2 +z ( 1 _ a ,4 )(1 _ x 4 ) ' 



n=0 

' ' 1 ' r/.f" | f/.r jf/./-;> 



n=0 ^ 



A).A)./o v/(l - X^X 2 ) 2 + z(l - x\)(\ - x\x%) ' 
The special values of (q(s) at s — 2, 3 are given as follows. 

Theorem 3.2 (Ichinose-Wakayama [8]). If a/3 > 2 (ie. < 1/(1 - a/3) < 1), then 



2 / . -, 2 



C (2) = 2 - ° + f „ a2 ,i, + fc|Y 92 f 1 



2y/aP(a(l-l) J y y, 2> \a + [3 J »°\af3 - I 
where £(s, x) = X)^Lo( n + 2; )~ s * s Hurwitz zeta function. 
Remark 3.3. We can determine the functions 32(3?) and 33(x) as follows: 

92 (x) = J 2 (0)g 2 (x) , 33 (x) = J 3 (0)g 2 (x) + 33 (x) , 

where 

. , . :i / ; :i ,r v /1 3 

32 W := 





1 




a 




; i + x) = 









v^T^ V 4 ' 4 ' 'l + z/ V 4 ' 4 

9 J21 / !\ 3 / \ n n— 1 -. / 1\- 2 

«"^»-'r 7 (rib) SftTTi - ' : 



See [18] and [11] for detailed calculation. ^> 

4 Higher Apery-like numbers 

Looking at the definition of J2(n) and Js(n), it is natural to introduce the numbers Jfe(n) by 

Jfc(n) := 2 k [ 
J [0,1 



(1 — xf)(l — x| • • • x^) \™ dx\dx2 ■ ■ ■ dxk 

V (1-X?..-X2)2 J l-X?...^ • 



We refer to Jfc(n) as higher Apery-like numbers. In fact, the generating function 



9k(z) := V f J fc (n)*» = 2* / dXl ^ 2 i i i (4.1) 

and its further generalizations are used to describe the 'higher' special values Cq(^) > 4 ) like Theorem 3.2 
(see Remark 4.1 below). 

It is immediate that Jfe(0) = (2 fe — 1)C(&)- Further, as we mentioned in [11], the formula 

3 Lfe/2J_1 i / i\ 1 - 
■MD - I E ^(*-^l)+-^ L (4-2) 
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holds (see §6.2 for the calculation). It is directly verified that 

4J fc (l) - 3J fc (0) = J fc _ 2 (l) (fc>4) 

Remark 4.1. We can calculate that 

2y/aP(a/3-l) J ^ v ' ' ' ' V« + /V ""V"/ 3 - 1 
+ 2 



fa(4)=2| B/ a , + ; j fc(4,i/2 )+ 4f^y 1 

2 

a — p\ /" \&dxidx2dx^,dx^ 



+ 



+ PJ J[o,i] 4 ^(1 - xfx|x§x|) 2 + 71(1 - xfx|)(l - 
a — /3\ 4 /" \6dx1dx2dx3dxi 

^TPJ J m * V(l - ^PPfP + 7i(l - - 4x1) + 72 (1 - xf)(l - xt)(l - * 4 )(1 - x\) 



where 71 = l/(a/3 — 1) and 72 = af3/(af3 — l) 2 . {} 

Similar to the case of J2(n) and J^{n), the higher Apery-like numbers Jfc(n) also satisfy a three-term 
recurrence relation as follows. 

Theorem 4.2. The numbers Jfc(n) satisfy the recurrence relations 

An 2 J k {n) - (8n 2 - 8n + 3)J k (n - 1) + 4(n - l) 2 J k {n - 2) = J fc _ 2 (n) (4.3) 

/or n > 2 and fc > 4. 

We give the proof of Theorem 4.2 in §5. It is remarkable that the left-hand side of (4.3) has a common 
shape with those of (3.1) and (3.2), and (4.3) gives a 'vertical' relation among Jfc(n)'s, i.e. it connects Jfc(n)'s 
and J/c_ 2 (n)'s. 

Example 4.3. First several terms of J 4 (n) are given by 

MO) = 15<(4), J 4 (l) = M C (4) + A C (2), J 4 (2) - ^C(4) + ^C(2), 

T 2205 3745 T 129735 1044135 

J 4 (3 = C(4) H C 2 , J 4 4 = C(4) H C(2 ), • • • 

KJ 256 sv ; 4096^ v h y> 16384 sw 1048576 sw ' 



We also see that 



4J 4 (1) - 3J 4 (0) = ^C(2) = J 2 (l), 
123 

16J 4 (2) - 19J 4 (1) + 4J 4 (0) = — C(2) - J 2 (2), 

64 

36J 4 (3) - 51J 4 (2) + 16J 4 (1) = ^C(2) - J 2 (3), 

25947 

64J 4 (4) - 99J 4 (3) + 36J 4 (2) = y^C^) - J 2 (4). 







Define another kind of generating function for J k (n) by 

~ f 1 f 1 f 1 1 - x 2 ■ ■ -x 2 



Theorem 4.2 readily implies the 
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Corollary 4.4. The differential equation 

„ Wk-2(t) - W k - 2 (0) . . 

T> H w k {t) = (4.5) 

holds for k > 4. Here T> H is the differential operator given in (3.7). 
Put 

Jo(n) := 0, Ji(n) := JJ fe (n) := J fe (n) (fc > 2). 

By Theorem 4.2, we have 

8n 3 J fc (n) - (1 - 2n)(8n 2 - 8n + 3)J fe (n - 1) + 2(n - 1)(1 - 2n)(3 - 2n)JJ fe (n - 2) = 2nJ fe _ 2 (n) 
for fe > 2 and n > 1. Hence, if we put 

a? d 2 d 

V w := 8z 2 (l + zfj-^ + 24z(l + z){\ + 2z)-^ + 2(4 + 27z + 27z 2 )— + 3(1 + 2z), (4.6) 

then we have the following (See also [11, Proposition A.3]). 
Corollary 4.5. The differential equations 

T\g 2 (.z) = 0, 
T\gz{z) = - 



1 + z' 



V w9k (z) = 2zl( 9k -^- 9k -^ ) (fc>4) 
hold. □ 

5 Proof of Theorem 4.2 

5.1 Setting the stage 

Assume k>2. We notice that 

n°° r°° e -(ti-\ — l- t fc)/ 2 fl — e~ 2tl ) n (l - er 2 ^ — I-**))" 
-y. ti-.-...^.,)^. — * " ' ** 

oo — u/2 /■« j.fe— 2 

**/ o ^(l-e-^l-e-^r* 



/ (l- e -«)2«+i y (fc_ 2 )! 
for each n > 0. Let us introduce 



Ii% = I { n %(u) ■= I 7^(1 - e- 2 r(l e-***Tdt 

for n, to > 0. We also put 

In k m(u) is symmetric in n and to if fc = 2 so that ln}m{u) = 0, but I^m(w) ^ in general. 
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It is convenient to set li%(u) = when k < 2. We sec that 

Thus we also set Jfc(n) = if k < 2. Under these convention, the following discussion for Jfc(n) is reduced to 
the one given by Ichinose and Wakayama [8] when k = 2, 3. 
For later use, we define 

a£\u):=l£l(u)=lW(u) (n>0), 

== \ + = tl-i (« > 1), 



so that 



5.2 Recurrence formulas for In)n{u) 
Integration by parts implies 

^' > = r(s(S)r) (1 - e "' r(1 - e "* +2 ' r<!i 



when n,m> 1. Since 



'o*^ 1 - e - M )"(^ 1 - e "^ m ) A 



|.(l_ e -2^ = 2ne- 2 *(l- e - 2 *)™- 1 



= 2n((l-e- 2 *)"- 1 -(l-e- 24 )"), 
— (1 - e - 2 "+ 2 *) m = -2me~ 2u+2t (I - e - 2tl + 2t ) m - 1 

2m ((1 - e - 2u+2t ) m - 1 - (1 - e - 2 «+ 2 *)™) 



for n, m > 1, we obtain the 

Lemma 5.1. The following three relations hold: 

ll^^in-m^-nl^ + ml^ (n,m>l), (5.2) 

n/W, - (2n - 1)/^ M „ + (n - l)/^ - me" 2 "!^,^ = \ (i^ 1^%) (n > 2, m > 1), (5.3) 

mi^ - (2m - l)4 fc Li + (m - 1)/£L 2 - ne" 2 "!^!,™-! = ^ - 4^) (n > l,m > 2). (5.4) 

□ 
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Plugging (5.2) into (5.3), we get 

li% ~ (I { n\m + + (1 " ^4^-1 = (n > 1, TO > 1), (5.5) 

which is a generalization of (4.14) in [8]. In particular, if we let n — to in (5.5), then we have 

4$ - + (1 - e- 2 ")/^,^ = 0. (5.6) 

Letting to = ra — 1 (or n = to — 1 and exchanging to by n) in (5.5), we also have another specialization 

" (if-Vn-i + C- 2 ) + (1 - e- 2u )li k \ n -2 = (n > 2), 

- A„ + 4-Vn-l) + (1 - e- 2U )I { n\n-l = (« > 2). 

Adding these equations, we get 

li fc L 2 = - &(«) + (1 - e - 2u )^!(«) (« > 2). (5.7) 

By specializing to = n in (5.3) and (5.4), we have 

n ^ri)n ~ (2" — 1)4-1,71 + ( n — 1)4-2,71 ~ ne 2 "4-l,n-l = l^{^n,n ^ — 4-l,7l)> (^-8) 

n/W - (2n - 1)4^ + („ - 1)JW_ 2 - ne-^Ln-i = ^-i - 4^) (5-9) 
for n > 2. Similarly, specializing m = n — 1 in (5.3) and n = to — 1 in (5.4) (and exchanging to by n), we have 



n ^n,n-l (2 n l)-^n-l,n-l + ( n lKra-2,n-l ( n l) e "^n-l,n-2 — 2 ^n-l,n-l)' 

1 

2 ( 



- (2n - l^-Vi + (» - l)^i,„- a " (» " l)^ 9 "^-! = ^(4-~5-i - 4-~fi) 



for n > 2. Adding each pair of relations, we obtain 

2nli k l - 2(2n - 1)1^ + 2(n - 1)I^_ 2 - 2ne- 2 "^_ ) li „_ 1 = ij^, (5.10) 
2nl^ - 2(2n - l^w + 2(n - 1)(1 - e~^\ n _ 2 =f n k ~!_\. (5.11) 

The formulas (5.6), (5.10) and (5.11) are rewritten as follows. 
Lemma 5.2. The equations 

ai k \u) + (1 - e-^ailfu) = 26«( u ), (5.12) 
«<#>(«) - (2n - 1)&W(«) + (n - 1)I^_2 - ™^ 2 "«^i(«) = ^^M, (5-13) 
n&«( u ) - (2n - ljaf^n) + (n - 1)(1 - e" 2 ")^^) = (5.14) 

fto/d. □ 

As a corollary, we also get 
Lemma 5.3. The equation 

na™(u) - (2n - 1)(1 + e^Ja^W + (n - 1)(1 - e" 2 ") 2 ^^) = 6^ fc_1) («) (5-15) 
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Proof. If we substitute (5.12), then we have 

b^iu) = 2nbi k Hu) 2(2n lja^u) + 2(n - 1)(1 - e" 2 ")^*) 
= n (a<*>(u) + (1 - e-*>)a<*Uuj) 2(2n lja^u) 
+ (n - 1)(1 - e- 2 ") (a^t*) + (1 - e^a^u)) 
= naW(n) - (2n - 1)(1 + ^(u) + (n - 1)(1 - e - 2 «) 2 a« 2 ( U ), 
which is the desired formula. 

Here we give one more useful relation. Using (5.2) twice, we see that 

= -n (-I { n\ n - (n - l)/^ + n/^,^) + n (l<2»-i - n^i.n-i + (n - 1)4' 
= n (26W( U ) - 2na<£ 1 («) + 2(n - 1)I^_ 2 ) • 

Thus we have 

a^ 2 \u) = 8n ( 6 W( U ) - na^u) + (» - 1)I^_ 2 ) • 
Combining (5.7), (5.16) and (5.14), we obtain 
Lemma 5.4. The equation 

a^- 2) {u)=Anb\t 1) (u) 

holds. 

In particular, the formula (5.15) is rewritten as 

nai k \u) (2n - 1)(1 + r 2 ")^,^) + (n - 1)(1 - e" 2 ") 2 ^^) = ^~ 2 \u). 

5.3 Relations for Bn\u) 
In view of (5.1), the differential 

^ag\u) = 2n f U J^-Al - e - 2t ) n e- 2u+2t {l - e~ 2u+2t ) n - x dt 
du Jo (fc-2)! 

is written in two ways as 

— oS n \u) = 2n (^In n-l ~ Inji) = ~2n (i^Ji ~ In-l,nj + ^n,n ^ 

for n > 1. Hence it follows that 

= -na( fe )( U ) + n(l - e- 2 >W(u) + la^u). 
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Using this formula, we have 



^Al k) (u) 2nsmhuA { n k l 1 (u) = e nu ^± a ( n k \u) + na ( n k \u) n(l - e" 2 ")^^)) 



= e ""^- 1 )( u ) = i-^- 2 )( u ). 



This implies the 

Lemma 5.7. TTie equation 



UK-D^W + i i- 



+ 2^ Jfe ~ 2 ^ -^fc-aW)- 



(5.20) 



Thus we obtain the 
Lemma 5.5. TTie equation 

2tanh||-B( fe )( U ) = SnB^u) - (2n + 1)B(*>(«) + tanh | *)(«) (5.21) 

holds for n > 1 . □ 
Remark 5.6. The differential of aQ fc ^(u) is given by 

%?><«> = 



v ""' (fc-2)! 

when k>2. If fc > 3, this is equal to 1 ^(m). 
We also see from (5.18) that 

nA^{u) - 2{2n - 1) coshu^^u) + 4(n - 1) sinh 2 u^ fc _> 2 (u) 



(5.22) 



(5.23) 



holds for n>2. □ 

5.4 Recurrence formula for Jk(n) 
Define 

i roo K( k )( v ) 1 /-oo 

K k( n ) = ^TTj o ^T du - M ^ = ^TTj Q tanh-B^^du. (5.24) 
By integrating (5.21) and (5.23), we have 

Jffc(n) = (2n + 1) Jfc(n) - 2nJ fe (n - 1) - M k (n), (5.25) 
2n( J fe (n) - K k {n)) = (2n - 1)(2 J fc (n - 1) - X fc (n - 1)) - 2(n - 1) J fc (n - 2) 

1 . , , „ , „ (5-26) 



Plugging these equations, we obtain 
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Lemma 5.8. Put 

L k {n) := J k - 2 {n) - J k - 2 (n - 1) + 2nM k {n) - (2n - l)M k (n - 1) - ^M fe _ 2 (n). (5.27) 
The recurrence formula 

4n 2 J k (n) - (8n 2 - 8n + 3)J k (n - 1) + 4(n - l) 2 J k (n - 2) = L fc (n) (5.28) 
ZioZds /or k > 2 and n > 2. □ 

When k = 2, the inhomogeneous term £2 (ft) in (5.28) vanishes and we get (3.1). When k = 3, we see that 
L 3 (n) — 2nM 3 (n) — (2n — l)M 3 (n — 1), which is equal to t^zwt (Lemma 6.3 in [8]), so we have (3.2). 

5.5 Calculation of the inhomogeneous terms 

Let us put 

*<">=- 14*- (5 ' 29) 

This definite integral converges if fc > 3. 
From (5.21), we have 

It follows then 

= 8n • 2 2 "" 1 Q fe (n - 1) - (2n + l)2 2n+1 Q fe (n) + 2 2 " +1 J fe _i(n), 

and hence 

Jfc_i(n) = (2n + l)Qfc(n) - 2nQ fc (n - 1) (5.30) 

for fc > 3 and n > 1. 

From (5.22), we also see that 

ntanh H _ 2 (2n - 1) f — ^- + tanh B { n k \(u) + 16(n - = J_ tanh " B^" 2 )( u ). 

2 " v ; v ; Vtanh| 2 J n ~ n ' K ' tanh § 4n 2 ™ w 



Thus we have 



which implies 



n2 2n+1 M fc+ i(n) - 2(2n - 1)2 2 "- 1 (Q fc (n - 1) + M k+l (n - 1)) 

+ 16(n - l)2 2 "- 3 Q fc (n - 2) = -^-2 2n+1 M k ^(n), 

An 



2nM k+1 (n) - (2n - l)M k+1 (n - I) - -^-M k -x{n) 

2n (5.31) 



= (2n - l)Q fc (n - 1) - 2(n - l)Q k (n - 2) 

for fc > 3 and n > 2. 
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Using (5.30) and (5.31), we obtain 

2nM k (n) - (2n - l)M k (n - 1) - — M fe _ 2 (n) 
= (2n - l)Q fc _i(n - 1) - 2(n - l)Q fc _i(n - 2) = J fe _ 2 (n - 1) 

for k > 4 and n > 2. Hence the inhomogeneous term is computed as 

L k (n) = Jfc_ 2 (n) - Jfe_ 2 (n - 1) + Jfc_ 2 (n - 1) = Jfc- 2 (n) (5.32) 

for fe > 4 and n > 2. This completes the proof of Theorem 4.2. 
Remark 5.9. It may be "natural" to assume (or interpret) that 



J {n) = 0, Ji(n) = 2 / (1 - x 2 )™" 1 ^ = 

Jo 



2"(n- 1)! 
(2n- 1)!! 



and 



uto(t) - 0, " Wl (i) - Wl (0)" = f] ^-1)^ " = 2t 2 F 1 (l, 1; I * 
Under this convention, Theorem 4.2 and Corollary 4.4 would include the case where k = 2, 3. {} 

6 Infinite series expression 

We give an infinite series expression of Jk(n). Using it, we prove the equation (4.2). 

6.1 Infinite series expression of J^{n) 
Let us put 

Us,t) := — J™ ( (1 " a4)( 2 1 ~ f4) y= (1 - s4 )"(! - i4 )"(! - ^ 2 )" 2 "- 1 - 



(1 - s 2 t 2 ) V (1 - s 2 t 2 ) 2 

Then we have 



Since 



J k {n) = 2 k / / ••• / fn(xi 1 x 2 ---Xk)dx 1 ---dxk- 
Jo Jo Jo 

us, t) = (i - , 4 )»(i - * 4 r E (~ 2 ™ " ') (-^ 2 )< = ^ E a + i)*. * 2i a - s4 )" <2i ( 1 - * 4 )"> 



it follows that 



2k 00 

Jk (n) - T^ry E (* + 1 ) 2 « J i (*> n)h-i(l, n). 

^ ''1=0 



Here I p (l,n) is given by 

fl r l rl 



I p (l,n):= / / •••/ («i • • • u p ) 2? (l - (ui ■ ■ -Up) 4 )"dui • • ■ du p . 
Jo Jo Jo 
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Notice that 

4 n n! 



h{l,n) = [ u 2l (l-u 4 ) n du 
Jo 

I p (l,n) = Y / (-iy( n ) / / ••• / (u 1 ---u p ) 2l+i idu 1 ---du p = V 

• =0 \Jy Jo Jo Jo • 



(2Z + l)(2Z + 5)---(2Z + 4n+l)' 

1 



(-1) 



Thus we obtain the expression 



j J {2l + 4j + l)P 



j f ,_2 fc 4"n!^ (l + l) 2n ^fn\ (-1)' 

H j (2n)! ^ (2Z + l)(2; + 5)---(2l + 4n+l)^U/(2^ + 4j + l) fe - 1 ' 1 ' J 



6.2 Example: calculation of Jfc(l) 
When n = 1, we see that 

7 m = 2 ([+lK[+2) / 1 1 \ 

& (2Z + 1)(2Z + 5) V(2? + l) fe - 1 (2Z + 5) fc -V 

= ~ (/ + l)(Z + 2)((2* + 5)*- 1 -(2/ + l)*- 1 ) 

^ (2Z + l) fc (2Z + 5) fc 

Using the identity 

(I + 1)(Z + 2) ((2/ + S)*- 1 - (2Z + l)^ 1 ) 

fc-2 

= ((2/ + l)(2l + 5) - (2Z + 1) + (2/ + 5) - 1) ]T(2/ + iy (21 + 5) k ~ 2 - j , 

i=o 

we have 

fc-2 

J fc (l) = 2 • 2 fc ^{5(fc - j -l, 3 + l)- 5(fc - j - 1,3 + 2) + 5(fc - 3,3 + 1) - 5(fe - 3,3 + 2))}, 



where 



Since 



5(a,/3) :=]T(2Z + l)-«(2; + 5)-' 3 . 



fe— 1 OO fe— 1 

53 su, k 3) = E E( 2Z + i) -7 '^ + 5 ) JW£ 

j=l i=0 j=l 

= g 1 2Z + 5 1 - (i$f ) fe_1 



(2Z + 5) fc 2Z + 1 1 _ ^ 2i+5 

i= 1 ^ 2Z + 1 



we have 



Jk(l) = 2 k+1 ( 



lA/ 1 1 \ 1 + 

~ 4 £jj ^ (2Z + l)^- 1 ~ (2Z + 5)*- 1 



2 + S(k, 1) - 5(1, fe) + 5(fe + 1,1) + 5(1, fe + 1) ) . (6.2) 



3/c+i 
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Let us calculate S(k, 1) and 5(1, k). By the partial fraction expansion 

1 ./. IN £ 



we see that 



x(x + a) k a k \x x + al ^— ' a fe_m+1 (x + a) r 



1 /lX^/l 1\ 1 fe / m + 2 1 



4 ) V 2Z + 1 2Z + 5 ) + 2 fe E ( 2 



(2/ + l) fe (2Z + 5) V 4 / \2Z + 1 21 + 5J 2 k r j^ 2 \ 2J (i + 3)" 1 ' 

1 1 \ 1 * /j\ fe ~" l + 2 ^ 



4) \ 2/ + l 2Z + 5) + 2 fe E (2 



(21 + l)(2l + 5) k V 4 / V2/ + 1 2Z + 5 7 2 fe ^ 2 V 2 / (i + 2 + |) r 

Thus it follows that 

k / 1 \ fc— m+2 1 / 1 ^ ^ 



^ 1 )= 2 rE(-2j ^'2 )+ 3(-4 

m=2 v 7 v 

t Gr + v. i> + 1 (if + 1 + 

m=2 v 7 v 7 

If we substitute these to (6.2), then we have 

(r, -1 k / , \ fe-m+2 -. 

m=2 v 7 

m=2 v 7 



3 V 4 

1 1 

3 ~ 3^ 



1 ^ / l\ fe -" l+3 ^ 1 1 / l^ fe 
+ ^TTE( v -2j « m '2 ) + H"4 y 

m=2 v 7 v ' 

fc+1 / -1 \ fc— m+3 ^ /lX' 1 1 



m=2 v 7 v 

Now it is straightforward to see that 



+ 77 + 



W = 3E^fe^C(m,i) + l±^i ) 

771 = 2 

_3 ^ 2" - 1 l + C-l)"- 1 

~ 4 2 fe - m 2 fe ~ 1 

2<m<k 
2 I fe-m 



fe-1 



3 Lfe/2J-i 
4 

m=0 



E"2- 2m c(^2 TO ,i) + i^#. 

m— n V / 



7 Differential equations for generating functions 

Utilizing the differential equations for the generating functions itffc(i), we give another kind of relations among 
the higher Apery-like numbers Jfc(n). 
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7.1 Equivalent differential equations 

Consider the inhomogeneous (singly confluent) Heun differential equation 

V K w{t) = u{t) 

for a given function u(t). Put z — and v(z) = (1 — t)w(t). Then we have 

V v{z) = —^—u 



z-1 \z-l 

Here V is the hypergeometric differential operator given by 

«.-*0-">£ + <i-*>s4 

We also remark that this is also the Picard-Fuchs differential operator for the family y 2 = x(x — l)(x — z) of 
elliptic curves. 

7.2 Recurrence formula for Jk(n) 

Put z = 4t an( i v k(z) = (1 — t)wk(t). By Theorem 4.2, Ufc(z) satisfies the differential equation 

1 °° / 



=f:fjE(-i) j+i (-)^(i+i)) 

n=0 \ j=0 V7/ / 



(7.1) 



The polynomial functions 
satisfy the equation 

(X>oPn)(«) = « n - (7-3) 
Hence we can construct a local holomorphic solution to (7.1) as 



(7.4) 



Notice that the difference Vk(z) — v(z) satisfies the homogeneous differential equation 

(V o (v k -v))(z)=0. (7.5) 

Thus it follows that 

v k (z) - v(z) = C k v 2 (z), (7.6) 

where the constant Ck is determined by 
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and v(0) is given by 
Therefore we have 

= (2 k - l)C(fc) 2 Fi Q, 1; + - «(0) 2 fi Q, 1; • (7.9) 

Consequently, we obtain the 

Theorem 7.1. WTien k> A, the equation 

JM -p-Kv)'0 (* - mk) -twrwdTp-O-^ + *>) (7 - io) 

holds. □ 
Remark 7.2. If we formally put Ji(n) = 2 ^n-ry.\ m (7.10), then we have 

«»)=B-')'(7)C) ('«>>-*£<5w(7 



since 



., 7 (2j + l)!! 2* + l 

This is nothing but the explicit formula (3.9) for J^{n). ^> 
Example 7.3. Since 

(2z + l) 2 



g(-iyQ« + i) = 3 ((2 )((7) -("*)) -3cm(-*) (i 

we have 



(2z + 2) 2 







7.3 Ascent operation and normalized higher Apery-like numbers 

For a given sequence {J(n)} n >o, we associate a new sequence 

«">" {§^(~<f <7 - n » 

Notice that J(O) 11 = 0. It would be natural to extend J(nf = if n < 0. By the discussion in the previous 
subsection, we have the 
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Lemma 7.4. Let {J(n)} be a given sequence and { J{ny} the one defined by (7.11). Then the equation 

An 2 J{nf - (8n 2 - 8n + 3) J(n - l) tt + 4(n - l) 2 J(n - 2) tt = J(n) (7.12) 

holds for n > 1 . 

Let us introduce the rational sequences Jfe(n) by 



Ji(n) 



2 n (n - 1)! 



2n — Yyyr (n>l), J 2 (n) := J 2 (n)/J 2 (0) (n > 0), 
Jfc(n) := Jk-2{n) (fc > 3, n > 0). 

We see that 

■/2 fc (l) = ^, W(l) = ^- (7-13) 

It is immediate to verify the 
Proposition 7.5. 

L fe / 2 J-1 l-(-l)k 

■h(n)= E C(fc-2m,l/2)J 2m+2 (n) + L_J_j fc („). (7.14) 



2 

m=0 



□ 



After this fact, we call Jfc(n) the normalized (higher) Apery-like numbers. By definition, Jfc(n) for k > 2 are 
written in the form 



p=0 



i\ 2 , 
■ ' n 



J k (n)=J2(-V P [ ~ 2 ) (J^(p), ( 7 - 15 ) 



where 



^ i '-^ 2 p - (i/2Mi)f r 

^ ( 2. + i)3 V zy (3/2)? r 



S a (p) = l, 5,(p) = -227^(7) ="2E 



p-1 , , lx -2 i 



(2i + l) 2 V i 7 ^ Vj 

i=0 v / \ / J=0 \J 

Thus it is enough to investigate Sfc(p) to obtain an explicit expression for normalized Apery-like numbers. 
Lemma 7.6. 

S fe+2 ( P + 1) - S k+2 ( P ) = ww - j—^. (7.16) 

Proof. By definition, we have 

+ 1) - s k+2 ( P ) = (^ /2 ) 2 ±(-iy Q J fe (, + 1). (7.i7) 

The sum in the right hand side is calculated as 

E(-i)'(*U(? + i) 

j=0 



{ n Sk(q) w i « j ■ § iij u + 1 ) ^ + 1} - 
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By the elementary identity 



("J (pj = Kp e z> ), 



we get 



Thus it follows that 

"i /a )*w-tiD'* (,,+i) 

- (2P+1) I, p J l(2p+l)» (2j, + 2)^' 

Therefore we obtain 

S k {j> + I) Skip) 



S k+2 {p + I) - S k+2 {p) = 



(2p + 2f (2p+l) 2 

as we desired. □ 

As a corollary, we readily have the 
Lemma 7.7. 

□ 

Using this lemma repeatedly, we obtain the 
Proposition 7.8. For each r > 1, 



\iiH Mr 



S 2 r+ 2 (p) = X! ( -2 ' gii,-,ir. (7-19) 

l<U<---<V<2p l--- r 

, , i /_i /2\~ 2 f_iVi+---+v 

where 

£ . = |0 1 < 3j < r s.f. = = 1 (mod 2), 2 ^ 

lv "' r 1 1 otherwise. 

□ 
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Example 7.9. We have 



s 4 (p) = 2^— — - 

i=i J 

*6lP)- 2^ i2 ,-2 £ *J ~ Z^ ,2,2 +2^727v*' 
l<i<j<2p J l<i<j<2p J i=l v ; 

58 (?) = Z^ ,2 ?2fc 2 £ '.J.fc 

l<i<7<fc<2p J 

/ 1 \ / (-1)» 

Z^ ,2 -2 fc2 i-l (2i) 4 ^ k 2 

l<i<j<k<2p J \ l<2i<2p v ; / \ l<fc<2p 







Remark 7.10. We see that 



lim S 2 (p) = 1, lim S 4 (p) = - — , lim S 6 (p) = - — . (7.22) 



In general, we can prove that 



}»» 'W 2 (p) = -p_T- (7-23) 



p— *-oo 



See [13] for the proof. <> 

8 Congruence relations among Apery-like numbers 

In this section, we study the congruence relation among the normalized Apery-like numbers introduced in the 
previous section. 

8.1 Congruence relations for Apery-like numbers 

We give several congruence relations among Apery-like numbers. 

Proposition 8.1 ([11, Proposition 6.1]). Let p be a prime and n = no + nip+ ■ ■ - + rikp k be the p-ary expansion 
of n € Z> (0 < nj < p) . Then it holds that 

k 

■M") = JJ J2(nj) (mod p). 

□ 

The following claim is regarded as an analogue of Proposition 2.2. 
Proposition 8.2 ([11, Theorem 6.2]). For any odd prime p and positive integers m,r, the congruence relation 

J 2 (mp r ) = J2{mp r ~ ) (mod p r ), 
Mp r )p 3r = J3(p r - 1 )p 3{r - 1) (mod p r ). 

holds. □ 
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Proposition 8.3. For any odd prime p, the congruence relation 

p-i 

^J 2 (n) = (modp 2 ) (8.1) 

n=0 

holds. 

Proof. We see that 

p-l p-l n /o A 2 / \ /o-\ 2p ~ 1 / \ 

E*w = ES-')'*-' ff) (") -B-i>"»-'(?) Eft) 

n=0 n=0 j=0 VJ 7 V7/ j"=0 V J 7 n=j V7 7 



2jV 1 



i / V i / i + 1 



pE i6 - j \;Jjtt (modp2) 

since ( j') 2 is divisible by p 2 il < j < p. Notice that 



16"' ( f ) = 2 . J )l^ 2 ) (modp) 



for < j < p. Hence we have 



n=0 j=0 V - 7 / \ J / J ~ 

By putting n = 2^- and m = in the identity (see [6, Chapter 5.3]) 



fe>0 

we have 



^ , n + k\ fn\ (-l) k = mini / m\ 

k \kk + l + m K ' (m + n + iy\n)' { ' ' 



Hence we obtain the desired conclusion. □ 
Proposition 8.4. For each odd prime p, it holds that 

J 2 (E=L) = A 2 (Z±) (modp 2 ). (8.4) 
Here A 2 (n) is the Apery number for ((2). 
Proof. It is elementary to check that 

(,- r + , )s( _ 1)t( - f) j 1+p |-_L_j (nodp2) 
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for k = 0, 1, . . . , (p — l)/2. Using these equations, we easily see that both A 2 (^~) and J 2 {^-) are congruent 
to 

(P"l)/2 



EM)'(-f) 



EW=f-) (modp 3 ) (8.5) 



modulo p 2 . □ 
Remark 8.5. The following supercongruence 

A 2 (^)^X p (modp 2 ) 

holds if p is a prime larger than 3 (see [9]; see also [15, 29]). ^> 
8.2 Conjectures 

In the final position, we give several conjectures on congruence relations among normalized (higher) Apery-like 
numbers. 

Conjecture 8.6 (Remark 6.3 in [11]). For any odd prime p, the congruence relation 

n=0 V y 

holds. 

Remark 8.7. The conjecture above is quite similar to the Rodriguez- Villcgas-type congruence due to Mortenson 
[14] 

g ( 2 ;)V« = (^) (mod, 2 ). (8.6) 

We also remark that the following "very similar" congruence relation is obtained in [19]: 

b-i)/2 /0 X 2 „ (p-i)/2 , N / x 

E O"-^^ 1 ^"" £ 0)Ht) ,modp3) ' (8J) 

where p is an arbitrary odd prime number. <0> 
The following conjecture is regarded as a "true" analogue of Proposition 2.2: 

Conjecture 8.8 (Kimoto-Osburn [10]). For any odd prime p, the congruence relation 

J 2 {mp r - 1) = (-^j J 2 (mp r - 1 - 1) (mod p r ) (8.8) 

ftoZds /or any m, r > 1. 

Remark 8.9. When r = 1, the conjecture (8.8) is obtained by using the elementary formulas 

and Mortenson's result (8.6). 
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Conjecture 8.10. For any odd prime p and m,r£ Z >0 with m odd, it holds that 

j 2 (l_g___) _ \ p J 2 ( m P r 2 1 - 1 ) + (-l)^- 1 >/yj 2 ( m P r 2 2 - 1 ) = (mod^), (8.9) 
where A„ is given by 

oo oo 
n— 1 n— 1 

Further, the congruence (8.9) ZioWs modulo p 2r if p > 5. 

Notice that (8.4) is a special case of the conjecture above (see [15, 29]). It is remarkable that both A 2 ( ' np 2 ~ :L ) 
and J2( mp 2 ~ 1 ) satisfy exactly the same congruence relation ((2.3) and (8.9)), though they are not congruent 
modulo p r in general. 

Conjecture 8.11. For any odd prime p, the congruence relation 

p-i 

^J 2 fc(n) = -1 (mod/) (8.10) 

71=0 

holds for any k>2. 
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